Hudson's Theorem for finite-dimensional quantum systems 
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We show that, on a Hilbert space of odd dimension, the only pure states to possess a non-negative Wigner 
function are stabilizer states. The Clifford group is identified as the set of unitary operations which preserve 
positivity. The result can be seen as a discrete version of Hudson's Theorem. Hudson established that for 
continuous variable systems, the Wigner function of a pure state has no negative values if and only if the state is 
Gaussian. Turning to mixed states, it might be surmised that only convex combinations of stabilizer states give 
rise to non-negative Wigner distributions. We refute this conjecture by means of a counter-example. Further, we 
give an axiomatic characterization which completely fixes the definition of the Wigner function and compare 
two approaches to stabilizer states for Hilbert spaces of prime-power dimensions. In the course of the discussion, 
we derive explicit formulas for the number of stabilizer codes defined on such systems. 



I. INTRODUCTION 



A. General Introduction 



The Wigner distribution establishes a correspondence be- 
tween quantum mechanical states and real pseudo-probability 
distributions on phase space. 'Pseudo' refers to the fact that, 
while the Wigner function resembles many of the properties of 
probability distributions, it can take on negative values. This 
phenomenon has been linked to non-classical features of such 
quantum states (see Ref. 1 1] for an exposition of literature on 
that problem). It is naturally of interest to characterize those 
quantum states that are classical in the sense of giving rise to 
non-negative phase space distributions. 

For the case of pure states described by vectors in 7i = 
L^(]R), the resolution of this problem was given by Hudson 
in Ref. f^. Later, Soto and Claverie generalized Hudson's 
result to states of multi-particle systems (Ref. [3]). 

Theorem 1. (Hudson, Soto, Claverie) Let ip e L^(R") be a 
state vector. The Wigner function oftp is non-negative if and 
only ifip is a Gaussian state. 

By definition, a vector is Gaussian if and only if it is of the 
form 

where q,x £ IR" and 9 is a symmetric matrix with entries in 

It is our objective to prove that the situation for discrete 
quantum systems is very similar, at least when the dimension 
of the Hilbert space is odd. Before stating the result, we pause 
for a brief overview of its main ingredients: discrete Wigner 
functions and stabilizer states. 

The Wigner function ||4J] of a pure state ip G L'^ (R) is com- 
puted as 



w^ip, q) 
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^(9 -7^0^(9 + 7^0- (1) 



Equivalently, W^p is the (symplectic) Fourier transform of the 
characteristic function S^, which in turn is defined by 

'^^(P,q) = tr(u;(p,9)T|-!/;)(V;|). 

Here, q) = e^iv^-iP) are the well-known Weyl or dis- 
placement operators filil]. Partly triggered by the advent of 
quantum information theory, considerable work has been un- 
dertaken to explore Wigner functions for finite-dimensional 
quantum systems Two 
approaches might be identified in the literature on that sub- 
ject. The first one aims to cast the definition of the Wigner 
function into a form that can be interpreted for both continu- 
ous variable and discrete systems [9, 10, 11, 15J. The second 
approach - introduced by Gibbons, Hoffman, and Wootters 
in Ref. - focuses on the properties of Eq. ([T]i. The au- 
thors imposed a set of axioms which a candidate definition of 
a discrete Wigner function would have to fulfill in order to 
resemble the well-known continuous counterpart. 
We will argue that, for odd dimensions d. 



W^{p,q)^d-^ J2 



^{q~2-'Oij{q + 2-'0 
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is the most sensible analogue of Eq. judged in terms 
of either of these approaches. Here, p, q are elements of 
Zrf = {0, . . . and2"i = (d+l)/2 is the multiplicative 

inverse of 2 modulo d. Indeed, the definition given above is 
the discrete symplectic Fourier transform of the discrete char- 
acteristic function and will be shown to be the unique choice 
to mimic certain desirable properties of the continuous Wig- 
ner function. 

Stabilizer states were originally defined by Gottesman in 
Ref. 117] as the joint eigenvectors of certain sets of elements 
of the qubit Pauli group. Exceeding the case of qubits, sta- 
bilizer states for higher-dimensional quantum systems have 
been treated in the literature (see, e.g. Refs. 1 1^ 
Such states find manifold applications in quantum information 
theory, ranging from quantum error correction ll22ll to Cluster 
state quantum computation fl^. Although displaying com- 
plex features such as multi-particle entanglement |24], stabi- 
lizer states allow for an efficient classical description. In par- 
ticular, a quantum computer that operates only with stabilizer 
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states can offer no principal advantage over classical methods 
of computing 1 22i1 . The latter statement is sometimes called 
Gottesman-Knill Theorem. 

Using that language, we intend to show: 

Theorem 2. (Discrete Hudson's Theorem) Let d be odd and 

ip e L^{Wi2) be a state vector. The Wigner function ofip is 
non-negative if and only ifij) is a stabiUzer state. 

Given that ipil) 7^ for all q, a vector ip is a stabilizer 
state if and only if it is of the form 

where q,x ^ TL"^ and 9 is a symmetric matrix with entries in 

Theorem 12] should convey two central messages. Firstly, if 
the right definitions are employed, the continuous and the dis- 
crete case behave very similarly (even though the methods of 
proof are completely different). Secondly, it adds further evi- 
dence to what might be called a piece of folk knowledge in the 
field of quantum information theory: namely that stabilizer 
states are the natural finite-dimensional analogue of Gaussian 
states. 

The paper is organized as follows. We survey previous 
work on the subject in Section ITB] Section is devoted to a 
superficial, yet self-contained introduction to Weyl operators, 
characteristic functions, Wigner distributions and stabilizer 
states. The main theorem is proven in Section |III] Sections 
rvl to I VII] address various related topics. The results of these 
last three sections do not rely on each other. Concretely, we 
comment on the relation between stabilizer states and Gaus- 
sian states in Section ITVl we consider mixed states with posi- 
tive Wigner functions in Section IV] and use Section [Vll] for a 
discussion of Hilbert spaces whose dimension is the power of 
a prime. 

Readers interested only in the structure of the proof, but 
not in its full generality, are deferred to Ref. ll25ll . where a 
particularly simple special case of the main result is laid out. 

B. Previous Results 

Recently, Galvao et. al. took a first step into the direction 
of classifying the quantum states with positive Wigner func- 
tion 1I27I1 . To explain the relationship of their results to the 
present paper, we have to comment on an axiomatic approach 
to discrete Wigner functions and, further, on stabilizer states 
in dimensions that are the power of a prime number 

In Ref. 1 16], Gibbons, Hoffmann, and Wootters listed a set 
of requirements which should be met by any definition of a 
discrete Wigner function W . Denoting the dimension of the 
Hilbert space by d, their axioms amount to 

1. (Phase space) is a linear mapping sending operators 
to functions on a d x d lattice, called the phase space. 

2. (Translational covariance) The Wigner function is co- 
variant under the action of the Weyl operators (in the 
sense of Theorem|7]i. 



3. (Marginal probabilities) There exists a function Q{X) 
that assigns a pure quantum state to every line A in 
phase space. If ip is state vector, then the sum of its 
Wigner function along A must be equal to the overlap 

Let us call functions that fall into this class generalized Wigner 
functions. This term is justified, as the characterization does 
not specify a unique solution: for a d-dimensional Hilbert 
space, there exist d'^^^ distinct generalized Wigner functions. 
Note also that the construction has been described only for the 
case where d = is the power of a prime, because only then 
the notion of a line in phase space has a well-defined meaning. 

We turn to the second remark, concerning stabilizer states. 
Consider a composite system, built of n d-level particles. We 
are free to conceive it as a single d" -dimensional object. The 
two points of view give rise to different definitions of sta- 
bilizer states, the 'single-particle' one being starkly reduced 
as compared to the multiple-particle one. In Section [VTI] we 
show that the set of single-particle stabilizer states is strictly 
contained in the set of multi-particle ones. Indeed, the ra- 
tio of the respective cardinalities of the two sets grows super- 
exponentially in 7i. As an example, the generalized Bell and 
GHZ states 

d-"/2 \l) ® d-"/2 J2 \i) ® \i) ® 

i i 

arguably the best-known multi-particle stabilizer states, do not 
belong to the respective single-particle sets. 

The result of Ref. t27il concerns quantum states in prime- 
power dimensions that are non-negative with respect to all 
possible definitions of generalized Wigner functions. These 
states are shown to be mixtures of single-particle stabilizer 
states, as described above. The authors aim to establish nec- 
essary requirements for quantum computational speedup. In- 
deed, if the Wigner function of a quantum computer is pos- 
itive at all times, then it operates only with stabilizer states 
and hence offers no advantage over classical computers, by 
the Gottesman-Knill Theorem. 

Thus for the case of non-qubit pure states. Theorem |2] im- 
plies the results of Ref. [27] and goes further in two essen- 
tial ways. Firstly, it suffices to look at a single definition 
of the Wigner function, as opposed to d"''' generalized 
ones. Secondly, quantum computation and the Gottesman- 
Knill Theorem are naturally set in the context of multiple par- 
ticles. Our definition assigns positive Wigner functions to all 
multiple-particle stabilizer states, while Ref. 1,271 effectively 
relies on the single-particle definition ll4lll . On the other hand, 
our main theorem does not address qubits or mixed states, 
which Galvao et. al. do. 



II. PHASE SPACE FORMALISM 

The term phase space formalism encompasses the ideas 
and tools in relation to the Weyl representation, to be defined 
shortly. We will give a concise introduction in this section. 
Many of the results presented can be found in the literature. 
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but some, e.g. the CHfford covariance of the Wigner function 
in non-prime dimensions, seem to be new. 



A. Weyl representation 

We start by considering a d-dimensional quantum sys- 
tem, d odd. In its Hilbert space H, we choose a basis 
{\Q), . . . ,\d — 1)}, labeled by elements of 'Ed- Henceforth, 
Z(j will be referred to as the configuration space and abbrevi- 
ated by Q. 

The pivotal objects in the phase space formalism are the 
Weyl operators (also known as the generalized Pauli opera- 
tors), as constructed below. Let xil) = e"r*9. The relations 



x{q)\x) = \x + q), z{p)\x) = x{px)\x) 



(2) 



define the shift and boost operators respectively. The Weyl 
operators are given by 



w{p, q) = x(-2 V?) 



(3) 



for p,q,t 6 Q. The specific choice of phases will prove use- 
ful later on [42] . The set of Weyl operators is closed under 
multiplication, up to phase factors. Direct computation shows 
that the composition law is given by 



w{p, q)w{p',q') 



(4) 



)w{p + p',q + q'). 



The square brackets denote the standard symplectic inner 
product on Z 



2. 



where 



J = 



1 

-1 



(5) 



(6) 



We write w{v) = w{vp,Vq) for elements v = {Vp,Vq) G Z^^. 
The space V :— Q x Q with inner product given by Eq. (|5]l 
will be called phase space in the sequel, owing to its analogy 
to the phase space known in classical mechanics. 

The preceding constructing generalizes naturally to multi- 
ple particles. Indeed, the configuration space of an n-particle 
system is given by Q = S^'. Multiplication between two 
elements p,q E Q is understood as the usual inner prod- 
uct pq — J^iPiq-i- The Hilbert space is again spanned by 
{\q)}qeQ and the Weyl operators are defined to be the tensor 
products 

w{p,q) = w{pi,...,pn,qi,---,qn) (7) 
= w{pi,qi) ® ■ ■ ■ (E)w{pn,qn)- 

Equations (|4]i, Q remain valid in the multiple-particle setting, 
if we substitute the matrix J by its multi-dimensional version 



J = 







nxn -D-nxn 
-D-nxn Onxn 



We end this section with some miscellaneous remarks. 

A state vector Itp) can be identified with a complex function 
on configuration space by setting ip{q) — (<z|V'}- We will use 
both representations interchangeably. 

The continuous Weyl operators w{p^q) = e*^^^^'^'^-', 
p, (? G R fulfill exactly the same composition law as stated 
in Eq. (|4|i, if x is set to x{q) ~ and the other symbols 
are interpreted in the obvious way. In fact, Eq. (|4| is then 
equivalent to the fundamental Weyl commutation relations |@]. 
Having this analogy in mind, p and q will sometimes be called 
momentum and position coordinates respectively. 

For future reference, note the two simple relations 



{w{p,q)i}}) {x) = x{ 

ti w{p,q) = d"- 6p^o6q^o 



2 ^pq+px)'ip{x 



q). 



(8) 
(9) 



It remains yet to justify the name Weyl representation. For 
V E V,t E TLd, define t) ~ x{t)w{v). Equation ^ takes 
on the form 

W{vi,ti)w{v2,t2) = w{vi + V2, ti +t2 + 2^ [wi , ^2] ) . 

The set V x K^, equipped with the above composition law is 
called the Heisenberg group H{7L''^), the Weyl matrices con- 
stituting a unitary representation of H{'E^i) |6]. This point 
of view on Weyl operators will be needed only in Appendix 

MM 



B. Clifford group 

The Clifford group is the subset of the unitary operators 
that map Weyl operators to multiples of Weyl operators under 
conjugation: 



Uw{v)U'^ = c{v)w{S{v)) 



(10) 



for some maps c : V ^ € wd S : V ^ V UJ]. The 
structure of the Clifford group is described in the following 
theorem |43]. 

Before stating the theorem, we have to comment on a re- 
appearing issue: namely that things are more involved if d is 
not a prime number For prime values of d, 7Ld has the struc- 
ture of 3l finite algebraic field, TH^ is definite vector space and 
most of the intuitions we have about vector spaces continue 
to be true. Among the more severe deficiencies of the gen- 
eral case is the fact that not every element a of TLd possesses 
a multiplicative inverse modulo d. But even if the analogue of 
a theorem about vector spaces holds for non-prime values of 
d, it is often difficult to find a proof in the literature. Appen- 
dices IIX C] and HXD] contain a collection of statements of this 
kind. Less technically inclined readers will not loose much by 
skipping these sections. 

For the sake of clarity of language, we call functions / on Q 
which fulfill /(Aa + 6) = A/ (a) + /(&) linear, disregarding 
the fact that Q might fail to be a linear space. Similarly, a 
subset 5 of Q that is closed under addition and multiplication 
by elements of TL^ is referred to as a subspace. We define a 
function S to be symplectic if it is linear and preserves the 
symplectic form: [S* • , • ] = [ • , • ]. 
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Theorem 3. (Structure of the Clifford group) 

1. For any symplectic S, there is a unitary operator fJ-{S) 
such that 

/i(S') w{v) /i(S')^ = w{S v). 

2. ^ is a projective representation of the symplectic group, 
that is 

ti{SMT) = e*V(^T) 
for some phase factor e"^. 

3. Up to a phase, any Clifford operation is of the form 

U = wia)fi{S) 
for a suitable a Cz V and symplectic S. 

The representation fi is called the Weil or metaplectic rep- 
resentation Sin. Theorem [3] is could be called a discrete 
version of the celebrated S'fone- von A^eMmann Theorem |6]. Its 
proof is not essential for understanding the further argument 
and has therefore been moved to Appendix lIX Al 

Note that a Clifford operation is connected to a vector a 
and a linear mapping S. This should remind us of a well- 
known structure on linear spaces: affine transformations. An 
affine mapping A is of the form A{b) = Sb + a where S 
is an invertible linear operator and a a vector. Let us call A 
symplectic if its linear part S is. 

We will frequently use the 'dot notation' to define functions 
of one parameter; for example writing S ■ + a for A. 

Lemma 4. (Clifford group and affine transformations) The 
mapping 

• + a I— > w{a)fi{S) 

is a projective representation of the group of symplectic affine 
transformations. 

Proof. All we need to do is to compare the composition law 
of the affine group 

{S ■ +a) o [T ■ +h) = S{T-+b) + a 
= ST-+{Sb + a) 

to the composition law of the representation 

w{a)fi{S) w{b)fi{T) = w{a) fi(S)w{b)fi{S)^ KS)fi{T) 
— w{a)w{Sb)fi{S)fj,{T) 
cx w{Sb + a)^{S T) 

which proves the assertion. □ 

The correspondence established by the last lemma will find 
a very tangible manifestation in Section IIIDI when we will 
see that the Clifford group induces affine transformations of 
the Wigner function. 



C. Fourier Transforms 

Let Q — S^' and / : Q ^ (D be a complex function on Q. 
The Fourier transform of / is 

= hp) = \Q\-'^' E (11) 

In the course of the main proof we will be confronted with 
Fourier transforms of functions which are defined only on a 
subspace of Q. If d is prime, then any subspace of Q = Z^' is 
of the form Z^' , for some n' < n, so no new situation arises. 

For non-prime dimensions, however, subspaces may not be 
as well-behaved. Consider as an example {0, 3, 6} C Zg. The 
set is closed under addition and multiplication, but can clearly 
not be written as . 

To cope with this problem, we will cast Eq. (fTTI) into a form 
that is well-defined for functions / on more general spaces. 
The construction is presented below. It can be found in any 
textbook on harmonic analysis (e.g. Ref. [29]). 

A character of Q is a function C • Q ^ C such that 
C(a + b) — C(a)C(^)- '^ny character of Q is of the form 
C('Z) = x{xq) for an appropriate x G Q (see Appendix lIX CI ). 
We can hence conceive the Fourier transformation defined in 
Eq. (fTTT i as a function of the characters of Q: 

/(c) = igr'/'Ec(9)/(9)- (12) 

We denote the set of characters of Q by Q* . With these no- 
tions, Eq. ( fT2b defines a function Q* (D. If, now, S is any 
subspace of Q and / a function on S, the Fourier transform 

is well-defined. 

For / : V <C, we define the symplectic Fourier transform 

as 

(.F5/)(a) = 1^1-1/2 Ex(k &])/(&)• (13) 

b£V 

Finally, take a note that the normaUzation in Eqs. ( fTTT i and 
(fT2] i has been chosen in such a way that Parzeval 's Theorem 

11/11 = 11/11 holds, where = 1/(9) p. 

D. Definition and properties of tiie Wigner function 

Employing Eq. (|9]l in conjunction with the composition law 
Eq. (IDi, one finds that the Weyl operators {w{p, q)} form an 
orthonormal basis in the space of operators on H with respect 
to the trace scalar product d^" tr( • ^ • ). The characteristic 
function Ep of an operator p is given by its expansion coeffi- 
cients with respect to the Weyl basis: 

Sp(e,x)=d-"tr(«;(e,x)tp). (14) 
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We mentioned in the introduction that the continuous Wig- 
ner function is the symplectic Fourier transform of the char- 
acteristic function JSl The two latter concepts have been 
defined for finite-dimensional systems in the preceding para- 
graphs. We can now state, in complete analogy to the contin- 
uous case: 

Definition 5. (Wigner function) Let d be odd, Q — TH]^ for 
some n. Let V, Ti be as usual and let p be a quantum state on 
H. 

The Wigner function Wp associated with p is the symplectic 
Fourier transformation of the characteristic function "Bp. 

An explicit calculation yields, for all a ^ V, 

{TsEp){a) = d-^" ^x([a,^'])trM5)tp) 
bev 

= d-tT{(^d-J2xiW,b])wib)^^p) 

=■ d-"tr(A(a)p), (15) 

where we have imphcitly defined the phase space point oper- 
ator A{a) fiel]. 

Theorem|6]lists a selection of properties of the Wigner func- 
tion. For a more thorough discussion, the reader is deferred to 
Refs. il[ll. 

Theorem 6. (Properties of the Wigner function) 

7. The phase space point operators have unit trace and 
form an orthonormal basis in the space of Hermitian 
operators on Ti.. Hence the Wigner function of an Her- 
mitian operator is real, and further, the overlap 

vev 

and normalization relations 

Y,Wp{v)^tip 

V 

hold. 

2. For a pure state ip, the Wigner function :— 
equals 

3. When computing marginal probabilities, the Wigner 
function behaves like a classical probability distribu- 
tion: 

Y,W^ip,q) = \4,{q)\'. 



4. The multi-particle phase space point operators factor: 

n 

A{pi, . . . ,pn,qi, ■ ■ ■,qn) ^<^A^'\pi,qi) 

i 

(and hence so does the Wigner function). 

5. It holds that A{0)\q) = \ —q). In other words, the phase 
space point operator at the origin equals the parity op- 
erator. 

6. The Wigner function Wp a of an operator product is 
given by the ★-product { also known as the Groenewold 
or Moyal product ii^/).- 

Wp,(u) = {Wp^W„){u) 
:= d-^Y^Wp(u^v)W„(u + w)x{\vM)- 

v,w 

Proof. The proofs are all straight-forward; we give only hints 
on how to conduct them. It will be essential to recall the well- 
known relation 

J2 x{xy) ^ d'' Sy,o, (16) 

for all y G T^. 

Indeed, the first claim can be proven by using Eq. (fT6b 
together with the definition of the phase space point operators 
Eq. ( fTsl l. Employ Definition |5] and Eq. (fT6ll to establish the 
second assertion, which in turn implies the third one. Theorem 
16 141 makes use of the fact that xiPl) — Yii x{Pi1i)'i see also 
Section fvnl for a very similar and more explicit calculation. 
The validity of the fifth statement is best shown using Eqs. 
(ill,®. 

Let us lastly turn to Claim |6] We have noted that the phase 
space point operators form an orthonormal system. Hence we 
can expand an operator p in terms of its Wigner function as 
p — Wp{v)A{v). Substituting p and a by their respective 
expansions in Wpa{v) = d^"^ ir{A{v)pa) yields the desired 
formula with the help of Lemma |29] □ 

The following statement will be vital to the proof of the 
main theorem. It assigns an elegant geometric interpretation 
to the Clifford group. 

Tlieorem 7. (Chfford Covariance) Let U — w{a)p,{S) be 
a Clifford operation. Let p' :— U pU'^ for some Hermitian 
operator p. The Wigner function is covariant in the sense that 

Wp{v)^Wp,{Sv + a). 
Proof. We compute the action of the Clifford group on the 
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phase space point operators. 

w{a)n{S)A{b)fi{S)^w{a)'' 

vev 

= d-"Y,x{[b,v])w{a)w{Sv)w{a)'' 

V 
V 

= d-" E x([&,^-V])x([a,z;'])zi;(z;') 

v' -. — S V 

The claim follows by use of Eq. ( fTSl l. □ 

Our definition of the discrete Wigner function coincides 
with the ones used in Refs. 11511 . It is further equal to 

Leonhardt's version ['s'], up to a permutation of points in phase 
space; it corresponds to choice (a) in Ref. 1 12] and lastly to 
G = Z2 in Ref. UJ]. One can show that W, as defined 
here, fulfills the axioms of Ref. ilM 

which had been laid out 
in Section II Bl Put differently, it is an element of the set of 
generalized Wigner functions. Gibbons et. al. remarked in 
Ref. liT6ll that among the generalized Wigner functions, some 
stand out by their high degree of symmetry. In our language, 
this symmetry is an incarnation of the Clifford covariance es- 
tablished in Theorem |2] Naturally, it is now interesting to 
ask how much freedom is left in the definition of a Wigner 
function, once one requires Clifford covariance to hold. We 
show in Appendix IIX B I that the definition used here is virtu- 
ally unique in that regard. 

E. Stabilizer States 

Using the composition law of the Heisenberg group Eq. (01, 
it is easy to see that two Weyl operators w(ui),w(u2) com- 
mute if and only if [ui, V2] = 0. Now consider the image of 
an entire subspace M under the Weyl representation w. The 
set 

w{M) ^ {w{m)\m G M} 

consists of mutually commuting operators if and only if the 
symplectic form vanishes on M: 

[mi, 1112] — 0, for all G M. 

Spaces of that kind are called isotropic. Clearly, if M is 
isotropic, then the operators w{M) can be simultaneously di- 
agonalized. We will see that if |Af | = d", the eigenspaces be- 
come non-degenerate and can thus be used to single out state 
vectors in the Hilbert space. A subspace M of V is said to be 
maximally isotropic if its cardinality equals d". See Appendix 
IIX CI for a justification of that nomenclature. 



Lemma 8. (Stabilizer States) Let M be a maximally isotropic 
subspace of V. Let v ^ V. Up to a global phase, there is a 
unique state vector |M, w) that fulfills the eigenvalue equa- 
tions 

x{[v,rn\)w{rn) \M,v) = |M, w) 
for all m G M. 

Proof. Existence: It is elementary to check that 



\M\- 



E X{[v,m])w{m) 



(17) 



is a rank one projection operator fulfilling the eigenvalue 
equations. 

Uniqueness: According to Appendix IIX CI there are 
characters of A/, each giving rise to a distinct projection op- 
erator as defined in the last paragraph. Two distinct operators 
of that kind are mutually orthogonal, because they belong to 
different eigenvalues of at least one of the Weyl operators. But 
dim H = \Q\ = and thus there is no space for more than 
one-dimensional solutions to the given set of equations. □ 

The state vector \M, v) is called the stabilizer state associ- 
ated to M and v. For obvious reasons, one refers to the set 
of operators {x{[v,m]) w{m)\m G M} as the stabilizer of 
\AI, v). Due to the isotropicity of M, the stabilizer is closed 
under multiplication and thus constitutes a group. Occasion- 
ally, we write \M) for \M, 0). To specify a stabilizer state, we 
need to specify a maximally isotropic space M. This is best 
done by giving a basis {mi, ... , rrik } of M. It is convenient to 
assemble the basis vectors as the columns of a 2n x fc-matrix, 
which is generally referred to as the generator matrix. As the 
choice of a basis is non-unique, so is the form of the generator 
matrix. 

A stabilizer state \M) is a graph state if it possesses a gen- 
erator matrix of the form 



(18) 



where is a symmetric n x n-matrix Il24tl . The designation 
stems from the fact that d can be interpreted as the adjacency 
matrix of a graph. Many properties of \M) are describable 
in terms of that graph alone 02411 . Some authors require the 
diagonal elements ; to vanish (equivalently, no vertex of the 
graph should be linked to itself), but we will not impose that 
restriction. Note that there exist considerably more general 
definitions of graph states ifTgll . 

Obviously, we will be concerned with Wigner functions of 
stabilizer states. Lemma|9]clarifies their structure. 

Lemma 9. (Wigner functions of stabilizer states) The Wigner 
function of a stabilizer state \M, v) is the indicator function 
on M -\- V. More precisely. 



d" 



a G Af + w 
else. 
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Proof. The representation given in Eq. ( fTTI l of v) deter- 
mines the characteristic function 

S|M,„)(6) = d-"x(b,6])5M(6). 
We compute the symplectic Fourier transformation: 

{Ts^M,v)){a) = d'2"^x([a,6])x([«,&])^A/(&) 

b&V 

fceM 

= d^" (5M-L(a - «)• 

Where 

= e y| [m, w] = for all m S M} 

is the symplectic complement of M in y. But 7\/ is a max- 
imally isotropic space and hence M = (see Appendix 
llXCT l. □ 

In particular we know now that the Wigner function of sta- 
bilizer states is non-negative. The next sections are devoted to 
the proof of the converse. 



III. DISCRETE HUDSON'S THEOREM 

A. Bochner's Theorem 

Define the self correlation function 

K^{q, x) = ij{q + 2-^x)ij{q - 
and note that the Wigner function fulfills 

^^P^l) = (19) 

xeQ 

Fix a go G Q- Designating the function p i— > VK(p, qo) by 
W( • , go), Eq. ( fT9] l says that • , go) is the Fourier trans- 
form of if (go, • )■ Therefore, W is non-negative if and only 
if the d" functions K{qQ, ■ ) have non-negative Fourier trans- 
forms. 

In harmonic analysis, the set of functions with non-negative 
Fourier transforms is characterized via a theorem due to 
Bochner It is usually proven either in the context of Fourier 
analysis on the real line or else, in full generality, for har- 
monic analysis on - not necessarily abelian - locally compact 
groups. While the former statement is not general enough for 
our purpose, the latter is not easily accessible. However, it 
turns out that in the discrete abelian setting an elementary 
proof can be given. It is stated in the next theorem, along 
with a variation for subsequent use. 

Theorem 10. (Variations of Bochner's Theorem) Let M be a 

subspace ofQ. Let / ; A/ — > (D. /f holds that 



L The Fourier transform of f is non-negative if and only 
if the matrix 

A\^f{x-q) {x,qeM) 

is positive semi-definite. 

2. The Fourier transform of f has constant modulus (i.e. 
|/(x)| = const j if and only if f is orthogonal to its 
translations: 

(/,i(g)/) - /»/(a^-'?) =0 

for all non-zero q G M. 

Proof. The following computation is a variant of a well- 
known fact concerning circulant matrices. We claim that 
any character C of M is an eigenvector of A with eigenvalue 
A = |Af|-i/2/(C)- Indeed, plug ging in the definitions yields 

(AC)(x) = 5]A%C(9) 

q 

= J2fix-q)aq) 

= E/('?)C(9)C(^) 
= vw\f ioa^)- 

There exist characters and thus equally many eigenvec- 
tors of A. Therefore, A can diagonalized. All its eigenvalues 
are non-negative if and only if / is non-negative. 

By the same argument, A is proportional to a unitary matrix 
if and only if |/(g)| is constant. But a matrix is unitary if and 
only if its rows form an ortho-normal set of vectors. □ 

From here, the proof proceeds in two steps. Section UlIBI 
harvests Theorem 1 101 II to gain information on the pointwise 
modulus \ip{q)\ of a vector with non-negative Wigner func- 
tion. Building on these finding, we will analyze the properties 
of such Wigner functions in Section ITlI CI 

B. Supports and Moduli 

Lemma 11. (Modulus Inequality) Let -ip be a state vector with 
non-negative Wigner function. 
It holds that 

> |^(g-a;)||V'(g + x)| 

for all g, x G Q. 

Proof. Fix a g G Q. As is non-negative, so is the Fourier 
transform of K^(q, ■ ). Bochner's Theorem implies that the 
matrix A^y = K{x — y,q) is positive semi-definite which in 
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turn implies that all principal sub-matrices are psd. In partic- 
ular the determinant of the 2x2 principal sub-matrix 

/ K^{q,0) K.^{q,2x)\ 
\K^{q,-2x) K^{q,0) J 

[ _ 1^(9)1' ijiq + x)^{q~x) \ 

\^p{q + x)tlj{q-x) li^iq)? J 

must be non-negative. But this means 

mq)\'-mq + x)i>{q-x)\^>0, 

which proves the theorem. □ 

We will call the set supp i/' of points where a state-vector is 
non-zero its support. S = supp ip has the property to contain 
the midpoint of any two of its elements. Indeed, if a, 6 S 
S, then setting q = 2^^{a + b) and x — 2^^(a — b) in the 
Modulus Inequality shows that 

m2-\a + b))\ > |^(a)||V(&)| >0, 

hence 2^^{a + b) G S. Let us refer to sets possessing this 
quality as being balanced. 

The following lemma clarifies the structure of balanced 
sets. Recall that a subset A of F is affine if A = M + v 
for a subspace M and some vector v. An affine space is a 
subspace if and only if it contains the origin 0. 

Lemma 12. (Balanced sets) A subset S of Q is balanced if 
and only if S is an affine space. 

Proof We show the 'only if part, the other one being simple. 

As both the characterizations of balancedness and affinity 
are invariant under translation, there is no loss of generality in 
assuming that G S*. We have to establish that S is closed 
under both addition and scalar multiplication. 

Let a G S. We claim that 

2-^XaeS (20) 

for ain G IN and A < 2'. The proof is by induction on I. 
Suppose Eq. ( |20] | holds for some I. If \ < 2'+^ is even, then 

2~^-^\a = 2-\\/2)b G S. Else, 

2-'-iAa = 2-i(2-'^a + 2-'^a) g5, 

which shows the validity of Eq. (l20l i. 

There exists an integer I > d such that 2' = 1 mod d. 
Indeed, by Euler's Theorem, 2'^^'^^ — 1 mod d, where is 
Euler's totient function. So I ~ d(p{d) satisfies the require- 
ments. Inserting I into Eq. ( l20l i. we conclude that A a G 5* for 
all A < 2''. Thus certainly Aa G 5 for all X £ 'Ed and we 
have proved closure under scalar multiplication. 

If a, 6 G S then, by the last paragraph 2a, 2b G S and 
hence 2^^ (2a + 26) G S, establishing closure of S under 
addition. □ 

Lemma 13. (Constant Modulus) Let be a state vector with 
non-negative Wigner function. Then \il>{ -) \ is constant on the 
support ofip. 



Proof. Pick two points x,q G suppV" and suppose \ilj{q)\ > 

Letting z = x—g, the assumption reads |'0('z)| > IV'(9+^)I- 
The Modulus Inequality, centered atq + z, gives 

mq + z)\^ > |^(g)||V(q + 2z)|. (21) 

As supp'i/' is affine, we know that ^{q + kz) ^ for all 
k G Zd. Hence Eq. (ISTT i. together with the assumption implies 

\^{q + z)\^ > \^{q + z)\\ij{q + 2z)\ 
<^ \i:{q + z)\ > \ip{q + 2z)\. 

By inducting on this scheme, we arrive at 

|V'(<Z)I > Hiq + z)\ > \^iq + 2z)\ > ■•• 

and therefore \ip(q)\ > |V'('Z + dz)\ = \ijj{q)\, which is a 
contradiction. 

Thus \ip{q)\ < \ip{x)\. Swapping the roles of x and q 
proves that equality must hold. □ 

At this point, we have full knowledge of the pointwise 
modulus of a state vector with non-negative Wigner function. 
The phases of • ) are, however, completely unknown. The 
section to come addresses this problem indirectly, by studying 
non-negative Wigner functions. 

C. Non-negative Wigner functions 

To motivate the following, assume for a moment that ip has 
a non-negative Wigner function and further, that ilj{q) ^ for 
all q. Choose a qo E Q and consider the function iy( • , (7o)■ 
Lemma[T3]implies that K^,{qo, ■ ) has constant modulus and 
hence - by Theorem 1 10121 - W{- ,qo) must be orthogonal to 
its translations. Clearly, a non-negative function possesses this 
property if and only if it is supported on at most a single point. 

There hence exists a po G Q such that W{p, qo) cx Sp,pg. 
This observation starkly reduces the possible forms of positive 
Wigner functions; it will be generalized to state vectors with 
arbitrary support in the next lemma. 

Lemma 14. Let ^ be a state vector If is non-negative, 
then it is of the form 

W^{v) ^ d-'^ St{v) 

where T C V is a set of cardinality d". 

What is more, ifO^T, then the set of elements of T with 
vanishing position coordinates 

{(p,0)GTbGQ} 

is a subspace of V. 

Proof. Let S = supp V'- Again, we may assume that 5 is a 
subspace of Q, for else we replace ip by 'w{—s)ip for some 
s G 5. It follows that supp K^, ^ S x S. Indeed, 

K^{q,x) ^ ^ q±2-^x G S 

4^ q e S Ax e S. 
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Denote hy S-^ = {q e Q\sq = for all s € S} the orthog- 
onal complement of 5* 1.44,1 . We will adopt the common nota- 
tion [p] = p+S-^ for cosets of S^. It should be clear that [p] is 
nothing other but the affine space with directional vector space 
given by and base vector p. The set S* of characters of S 
can be identified with Q/ S^. Certainly, s ^ x{ps) defines a 
character of S for every p e Q. Further, x(ps) = x{p's) for 
all s G S* if and only if p — p' G S-^. That indeed all elements 
of S* can be obtained this way is shown in Corollary |26l 

Define K'^^ to be the restriction of to its support 5x5. 
For the rest of the proof, we fix a go G 5*. Now consider 

Wip,qo) = (i-"^x(pa;)if(go,x) 

xes 

Viewed as a function in p, W{p, go) has constant values on 
cosets of 5^. Therefore, 

W'{[p],q):=r\S\-^/^W{p,q) (22) 

is a well-defined function on 5*. The considerations of the 
previous paragraph allow us to identify W' {[■],qo) as the 
Fourier transform of K'{qQ,-). 

We can now repeat the argumentation presented just before 
the current lemma. Indeed, the modulus of K' {qQ,[-]) is con- 
stant and W is non-negative. Furthermore, by definition of 
qo, K'{qo,[-]) is non-zero and we may thus conclude that 
p 1— > qo) is supported on exactly one coset [po]. 

Normalization of implies ]'(/'( •)! = |5|~^/^. Hence 

X 

By Parzeval's Theorem, | |W^'([ • ], go)| P = l^l"^ as well. It 
follows that VK'([po],9o) = |5|"i/2. 
Inverting Eq. (l22T i gives 

Wip.,) = -'-'{l (23) 

which proves the first claim of the lemma. The cardinality of 
T is fixed by the normalization of the Wigner function (Theo- 
rem |6l6]l. 

Now suppose W{0,0) = W^'([0],0) ^ 0. Clearly, then 
W{p, 0) is non-zero if and only if p e [0] 44> p £ 5^. The 
last assertion of the lemma follows, since 5^ is a subspace of 

g. □ 

So a non-negative Wigner function is the indicator func- 
tions of some set T. This finding is compatible with Lemma 
|9l which describes the structure of Wigner functions of stabi- 
lizer states. The next two lemmas verify that T has indeed all 
the properties of the sets that appear in Lemma|9] 

Lemma 15. Let tp be a state vector. IfW^, is of the form 

W^iv) = d--^ St{v), 
then T is an affine space. 



Proof. The proof proceeds similar to the one of Lemma [12] 
There is no loss of generality in assuming that G T. 

First, we show that T is closed under scalar multiplication. 
To this end, pick a point a € T. There exists a symplectic 
mapping 5 that sends a to a vector a' of the form (a^, 0) where 
a'p £ Q (see Appendix BXDb . The set T' = 5 T is the support 
of the Wigner function of ii{S)'4!. By the second assertion of 
Lemmafltl Aa' G 5T for every A G Z^. Hence 5"i(Aa') = 
Aa G T. 

Turning to closedness under addition, let a,b £ T. By the 
last paragraph, 2a, 2b G T. Arguing as before, note that the 
set T — 2a is the support of the Wigner function of w{—2a)4' 
and thus closed under multiplication. As 26 — 2a G T — 2a, 
we know that b — a £ T — 2a and hence b + a £ T. □ 

Lemma 16. Let tp be a state vector such that is of the 
form 

W^{v)^d-''dTiv). 
IfT is a subspace, then it is isotropic. 

Proof. The vector tp describes a pure state, hence W^-kW^ = 
(recall the Moyal product, introduced in Theorem|6ll. Let 
u £ T. Plugging in the definitions gives 

= d-" W^{u + v)W^{u + w)x{[v,w]) 

v,w£T 

Note that J^wf^rXib^^w]) < \T\ = with equality if and 
only if [z;, w] =0 for all w. Hence 

For the left-hand and the right-hand side to be equal, T must 
be isotropic. □ 

Therefore T, as defined above, is of the form T — M + v 
where M is an isotropic space of cardinality rf". But then, 
Wti, is the Wigner function of a stabilizer state, by Lemma |9] 
We have proven: 

Theorem 17. (Main Theorem) Let ip G L'^{'E^) be a state 
vector If the Wigner function of tp is non-negative, then tp is 
a stabilizer state. 



IV. DISCRETE GAUSSIANS 

It has long been realized that the coefficients of stabilizer 
state vectors are described by quadratic forms. However, the 
current literature either neglects the non-prime case (Refs. 
||l9, 20, 26]) or is less expHcit (Ref. |21]) than the following 
lemma in showing the tight relation between Gaussian states 
and stabilizer states. 
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We will concentrate on stabilizer states with full support. 
This constitutes only a modest restriction of generality. In- 
deed, let i/jhea general stabilizer state, let Q' := suppi/;. Let 
us for the sake of simplicity assume that d is prime and Q' is a 
subspace of Q. The restriction of the coordinate function ip{q) 
to Q' can be thought of as defining a vector ip' of a quantum 
state of an n' := dimQ' particle system. It is now possible 
to check that ip' is a stabilizer state. In this way any stabilizer 
state can be viewed as one with full support, possibly on a 
smaller system. We will, however, not take the time to make 
this construction precise nor will we rely on it in this paper 

Lemma 18. Let ^ be a state vector. The following statements 
are equivalent. 

1. "0 is a stabilizer state and ijjiq) ^ 0/or all q G Q. 

2. Up to the action of a Weyl operator, ip is a graph state. 

3. There exists a symmetric n x n-matrix 9 and an x G Q 
such that 

Proof. (li]=>|2|. By assumption = \M,v) for some max- 
imal isotropic space M and a vector v. We claim that there is 
no non-zero p e Q such that (p, 0) e M. 
For suppose there exists such a p. Then 

{q\w{pMM) = xi-pqMM). 

On the other hand, 

(q|u;(p,0)|Af} =x(b,b,0)])(<?|M), 

by the definition of |Af, v). Hence supp \ M) must be con- 
tained within a hyper-surface of Q specified by pq = const, 
which contradicts the assumption that supp?/' = Q- 

There are d" elements in A/. By the last paragraph, no two 
of them have the same position coordinates. As there exist 
only d" = IQI possible choices for the position coordinates, 
one can find for every q £ Q a p e Q such that (p, q) G 
M. Let ei, . . . , e„ denote the canonical basis of ZJJ. Choose 
mi, ... , m„ e M such that the position part of rrii equals 
ei. The span of {m^li^i „ has clearly cardinality d", so 
we have found a basis of A/. By construction, the generator 
matrix composed of these basis vectors has the form shown 
in Eq. (fTsT l with some n x n-matrix 9. It is not hard to see 
that M is isotropic if and only if 9 is symmetric, establishing 
that \M) is a graph state. Theorem|2]and Lemma|9]show that 
w{v)\M) = \M,v) = IV^). 

(|2] ^ Let M be an isotropic space which possesses a 
generator matrix of the form given in Eq. ( fTST l. Let = 
Ci) be the ith column of that matrix. We need to establish 
the existence of a symmetric matrix 9 and an a; G Q such that 

{q\M,v) V((7). 

Indeed, choose 

6* = 2-1^9, x^ = [v,m,]. 




FIG. 1: Wigner function of the antisymmetric vector \'4>-)- 

Using Eq. (HJ, one can then check by direct computation that 
■4> fulfills the defining eigenvalue equations 

x([w,mj])w(mj)V' = ^ 

and hence = | A/, v), by Lemma|8] 

([5] =^ |7). Reverting the previous proof shows that is a 
graph state. It has maximal support by definition. □ 

The claimed analogy between stabilizer states and Gaussian 
states is apparent when comparing statement 3 to Theorem[T] 

V. MIXED STATES 

It is natural to ask how the results obtained before gener- 
alize to mixed states. Certainly, mixtures of stabilizer states 
are non-negative on phase space and it might be surmised 
that all such quantum states are convex combinations of sta- 
bilizer ones. In the context of continuous variable systems, 
Brocker and Werner refuted an analogous conjecture by giv- 
ing a counter-example 13 ill . Again, the situation is similar in 
the finite setting, as will be shown now. 

As a consequence of Theorem I6I5I A{Q) can be decom- 
posed as A{Q) = P+ + P-, where P± denotes the projector 
onto the symmetric and antisymmetric state vectors respec- 
tively. Since P++P_ = 1, we have that F_ = 1/2(1- A(0)). 
Because we know the Wigner functions of both 1 {W{v) = 
d^") and of ^(0) {W{v) = (5„.o), we immediately obtain 

^^-(^)=2\'^-" else. (24) 

For a single three-dimensional quantum system there exists 
a unique antisymmetric state vector = 2^^/^(| + 1} — 
— 1)), hence P_ = ("0- 1- Figure |2] depicts the Wig- 
ner function of the state p, obtained by mixing the pure states 
w(— 1, O)|'0-), w(— 1, with equal weights. 

The Wigner function of a single-particle stabilizer state is a 
line in the two-dimensional phase space, according to Lemma 
|9] There are d(d + 1) such lines and hence equally many 
stabilizer states. Assume these states have been brought into 
some order and denote the associated projection operators by 

Pi, ... , Pd(d+i)- Let p ^ Yl'^'"'^^^^ XiPi be a convex decom- 
position of p in terms these operators. If there is a point v in 
phase space where Wp{v) = and Wp- (v) ^ 0, then clearly 



11 



-1 1 




-10 1 
Momentum 



FIG. 2: Wigner function of the equal mixture of the vectors 

\il}—),w{—l, 0)\il)-) and ■w{ — l, — White squares stand for 
a value of 1/6, blaclc squares for 0. 

FIG. 3: The white squares mark all lines in S3 that do not intersect 
any point where the Wigner function shown in Fig. [2]vanishes. 

Xi must vanish. By exhaustively listing all 12 lines in Z§, one 
finds that p can have non-zero coefficients only with respect 
to the stabilizer states whose Wigner functions are shown in 
Figure |3] 

But p admits no convex decomposition in terms of these 
three lines. Indeed, no two of them cover all the points in the 
support of Wp, so only a mixture of all three lines could poten- 
tially suffice. Now notice that the point (1, —1) is an element 
only of the third line, while (1, 0) is contained in both the sec- 
ond and the third one. Therefore any mixture of these three 
lines takes on a higher value on (1, 0) than on (1, —1). The 
distribution Wp, on the other hand, is constant on its support. 

VI. DYNAMICS 

Having established which quantum states give rise to non- 
negative phase space distributions, the next step is to charac- 
terize the set of operations that preserve this property. We 
have seen in Section HID I that Clifford unitaries implement 
permutations in phase space and thus manifestly preserve pos- 
itivity. They are unique in that regard, as will be shown now. 

By the results of Section |III1 it is apparent that a unitary 
operation U can preserve positivity only if it sends stabilizer 
states to stabilizer states. One can reasonably conjecture that 
only Clifford operations possess this feature and in the case 
of single-particles in prime-power dimensions, a proof of this 
fact has been given in Ref . [27J . The general case, however, 
poses surprising difficulties which have forced us to take a less 
direct route. 

Let us shortly pause to clarify our objectives. We aim to 
characterize the set of unitaries U that satisfy statements of 
the kind: Wup^ is non-negative whenever Wp is. We can re- 
quire the above statement to hold for any Hermitian operator 
p, or just whenever p is a quantum state. In the former case 
the restrictions on U are much stronger than in the latter one. 



Indeed, by considering the image of the phase space point op- 
erators A{a) under the action of U and making use of Lemma 
|29l it is straight-forward to prove that only Clifford operations 
can preserve positivity of the Wigner functions of general Her- 
mitian operators. The following theorem is slightly more am- 
bitious in considering only the action of U on quantum states. 

Theorem 19. (Only permutations preserve positivity). Let 
U be unitary. If, for all quantum states p with non-negative 
Wigner function, it holds that Wijpu^ is non-negative, then U 
is Clifford. 

Proof. Firstly, take a note that substituting 'quantum state' by 
'positive operator' in the above theorem, only amounts to a 
change of normalization and does not alter the statement. Set 

p{p) := TmnWp{v), 

v{p) minargVFp :— {v G F|Wp(u) = p-{v)}. 

Let p be such that p{p) < 0. We claim that p{p) = p{p'), 
where p' — UpW. In other words: U preserves minimal 
values. 

Indeed, there exists positive constants X12 such that 

Xipip') + Aad-" = 0. 

Hence a := Xip + A2I has a non-negative Wigner function. 
The assumption p{p') < p{p) yields 

WuaW («) = + Asd-" < 

for every v £ v{p'), which contradicts the defining property 
of U . Thus ijl{p') < p.{p)- Substituting U by U^^ shows that 
equality of fj,{p) and p{p') must hold. 
Now set 

p(a) := (1 - w(a)P_ u>(a)1' 

for all a € V. We have /i(/3(a)) — p.{p'{a)) = — landt/(p) = 
{a}. The crucial observation lies in the fact that i'{p') contains 
only a single point as well. So, U preserves the 'pointed' 
shape of Wp{a). To see why that is the case, suppose there is 
a flo such that |j/(p(ao)')| > 1- There are d^" operators p(a)' 
and equally many points in phase space, so there exists an oi 
such that ^(ao) and v{ai) intersect in at least one point v. 
Define cr = l/2(p(ao) + p(ai)). It holds that ^(cr) > -1/2, 
whereas Wcr'{v) = —1 which is a contradiction. There is 
hence a well-defined function S which sends a to the unique 
element of i/(p(a)'). 

Finally, let a be any density matrix. The idea is to mix a 
very weakly to p{a), so that the positions of the minima of the 
mixture are still determined by p{a). Indeed, there exists an 
e > such that 

i^{p{a) + ea) ~ {a} 

p{p{a)+ea) = -1 + eW^ia); 

i^ipiaY + ea') = {5(a)} 

pip{ay + ea') = -I + eW^{S{a)). 

Hence Wa'{Sa) — W^ia). We have established that U acts 
as a permutation in phase space and is therefore Clifford by 
Lemma|29] □ 
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VII. PRIME POWER DIMENSIONS 

Wigner functions for quantum systems with prime power 
dimensions have received particular attention in the literature 
(most prominently in Ref. tl^T ). Once again, this is due to the 
fact that a finite field of order d exists exactly when d is the 
power of a prime and that the field's well-behaved geometri- 
cal properties facilitate many constructions. The present sec- 
tion briefly addresses the relationship between three natural 
approaches to Wigner functions for such systems. We assume 
the reader is already familiar with the definition of Weyl oper- 
ators over Galois fields; a thorough introduction can be found 
inRefs. Hill. 

Let d = for some prime number p. There are three 
natural ways of associating a configuration space to H.. These 
are 

1 . an n-dimensional vector space over Zp, 

2. a one-dimensional module over or 

3. a one-dimensional vector space over the Galois field 
Fpn of order p". 

The first and the second of these points of view have mani- 
festly been covered in this paper So far we neglected case 3, 
because - as we will see - it can be completely reduced to the 
first one. 

Let us quickly gather some well-known facts on finite 
fields. If p is prime and n a positive integer, Fpn denotes the 
unique finite field of order d = p". The simplest case occurs 
for n = 1, when Wp ~ Kp. For n > 1, fields Fp^ are realized 
by extending Fp, which is then referred to as the base field. 
Extension fields contain the base field as a subset. The exten- 
sion field possesses the structure of an n-dimensional vector 
space over the base field. A set of elements of Fpn is a basis if 
it spans the entire field under addition and Fp-multiplication. 
After having chosen a basis {bi, . . . , 5„}, we can specify any 
element / = f ^bi by its expansion coefficients {/*}. The 
operation 

n-l 
fc=0 

takes on values in the base field and is Fp-linear Therefore, 

^Tr(/5) 

defines an Fp-bilinear form. For any basis {bi}, there exists a 
dual basis {6*} fulfilling the relation Tr(6*6j) — 5i,j (we do 
not use Einstein's summation convention). From now on, we 
assume that a basis 6, and a dual one 6' have been fixed. 

Repeating the construction put forward in Section HIl we 
introduce the Hilbert space H = L^(Fpn ), in other words, Ti 
is the span of { | <z} | g G Fp^ } . The choice of a basis induces a 
tensor structure on H via 

i i 



We obtain a character of Fp^ by setting Xp" if) = 
Xp(Tr/). Note that for n ~ 1, Xp^ = Xp- Expanding mo- 
mentum coordinates p — J2j Pj^^^ the character factors: 

Xipq) = XpiJ^P^q' Tr(fo,6^")) ^l[xp{P^ql■ 

Similarly, the shift and multiply operators factor with respect 
to this tensor structure: 

i j i 

i j i j 

- ^z^'Hp^)\xi, 

i 

where x^*'' and z^*) act on the ith p-dimensional subsystem. 
A straight-forward computation along the lines just presented 
shows that both the Weyl operators and the phase space point 
operators factor: 

w{p,q) = (^w^'\pi,q')=w{pi,...,pn,q^,...,q") 

i 

A{p,q) = (g)AW(p„g*) = ^bi,...,p„,g\...,9"). 

i 

The above result thus states that the Wigner function induced 
by the choice Q = Fpn coincides - up to re-labeling of the 
phase space points - with the one for Q = F^. In particular, 
both definitions give rise to the same set of states with a non- 
negative phase space distribution. 

For stabilizer states, however, the situation is not as easy, 
as will be discussed subsequently. The preceding discussion 
suggests defining a map i : Fp„ — > Fp" by 

1-^ (pi, . . . ,p„,g\ ... ,9") 

(see Refs. H Q). Let M be a maximal isotropic sub- 
space of Fp„. It is readily verified that i{M) C F^" is again 
isotropic and a subspace. Further, we have shown that the sets 
of Weyl operators w{M) and w{l{M)) coincide and hence so 
do the stabilizer states \M) and \i{M)). 

The converse is not true. does not necessarily map Fp" 
subspaces to those of Fp„ . More precisely, if AI C F^" is 
a subspace, then t~^(Af) can easily be proven to be closed 
under addition, but will in general fail to be closed under 
Fpn -scalar multiplication. This proves the remark made in 
the introduction, namely that the set of 'single-particle' (i.e. 
Fp„) stabilizer states is a true subset of corresponding 'multi- 
particle' set. The following subsection gives a quantitative 
account of the relation of the sets. 

A. Counting stabilizer codes 

We are going to count the number of stabilizer states of a 
system composed of n d-level particles. In fact, the compu- 
tation given below is slightly more general in that it gives the 
number of /c-dimensional stabilizer codes ITtIi . 
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Stabilizer codes are generalizations of stabilizer states. Re- 
call Eq. ( fTTI l. where we showed that summing Weyl operators 
w{m) over the elements to of a maximal isotropic subspace 
M of V yields a one-dimensional projection operator. It can 
be shown that if the requirement of maximality is dropped, 
the sum still evaluates to a projector The range of this opera- 
tor is the stabilizer code defined by M. The dimension to of 
M and the dimension k of the stabilizer code are related by 
k = d"-™. 

Theorem 20. (Number of isotropic subspaces) Let V be a 
2n- dimensional symplectic vector space over Wd, where d is 
the power of a prime. The number of m-dimensional isotropic 
subspaces of V is given by 



Iso(n, TO, d) 



Hid" 



1), 



where the square brackets denote the Gaussian coefficients 



m — 1 

n 

1=0 



d"- 



- 1 



Proof. The proof is inspired by a method employed in Ref. 
ll33ll to solve a related problem. We count the number of lin- 
early independent m-tuples consisting of mutual orthogonal 
vectors. Indeed, as the first vector vi we are free to choose 
any non-zero element of V. There are d^" — 1 such choices. 
The second vector must lie in the symplectic complement of 
the span of the first vector (wi)^. Hence, V2 can be chosen 
from a 2n — 1-dimensional vector space, the only restriction 
being that V2 ^ (vi). It follows that there exist d^"^^ — d^ 
possibilities for 112. Inducting on this scheme gives 



i=0 



•2n~ 



(25) 



such tuples. 

However, since two different tuples might correspond to the 
same isotropic space, Eq. dZSl l over-counted the subspaces. 
To take that fact into account, we must divide by the number 
of bases within an TO-dimensional space. Arguing in a similar 
fashion as before, we arrive at HZlo^ ('^™ ~ fo'" sought- 
for number (see also Ref. [133]). Division gives 



Iso(n, m,d) — 



n 



Expanding d^^" — 1 = (d" * — 1) (d" ' + 1) and using 
the definition of the Gaussian coefficients concludes the proof. 

□ 

Corollary 21. The number of d"^™' -dimensional stabilizer 
codes defined on n d-level systems is 



Stabs(n,TO, d) = d" 



n(d"-'+i). 



In particular, the number of stabilizer states is 

71 

Stabs(n,n,d) = d" J|(d' + 1). 

i=l 

Proof. We only need to justify the pre-factor d"*. The defin- 
ing Eq. SvH generates a projector onto a stabilizer code given 
an isotropic space M and a character x(bi ']) on M. If 
dimAf = TO, then there are |M| = d™ distinct such char- 
acters (see Appendix lIXCI ). □ 

We can now compare the number of stabilizer states for n 
particles of dimension d to the corresponding number for a 
single d"-dimensional system: 

Stabs(n,n,d) _ HLiK + 1) ^ J| (-^^ _^ 



Stabs(l, l,d" 



d" + 1 



i=l 



> d^^=^' = d5 ("'-"). 



This is the super-exponential scaling mentioned in the intro- 
duction. 
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IX. APPENDIX 

A. Discrete Stone-von Neumann Theorem 

This section generalizes well-known results for prime- 
power dimensions (see e.g. Ref. ll34ll and citations therein) 
to all odd d. The proof is based on some simple observations 
employing group representation theory. We state a preparing 
lemma beforehand. 

Lemma 22. The Weyl representation is irreducible. 
Proof. We compute 

Y: |tru;(a,t)p = d-(^"+i) ^ 1 tr«;(0, t)^ 



aev. 
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which establishes irreducibility by a well-known criterion 
from group representation theory (see any textbook on that 
topic, e.g. [35]). □ 

Proof, (of Theorem^ By the composition law Eq. ^ it is 
clear that w'{p,q,t) :— w{S{p,q),t) is a representation of 
the Heisenberg group which affords the same character (i.e. 
tr w{a, t) — trw'(a, t)). The preceding lemma yields that w 
and w' are equivalent and thus the existence of follows. 
Further, 

^Ji{S)^l{T)w{p,q)^i{T)^ p{S)^ = ^i{S)w{T{p,q))^l{s)^ 

= w{ST{p,q)) 

= ^i{ST)w{p,q)^l{STy. 

Because the Weyl matrices span the set of all operators, the 
last line fixes /i(S'T) modulo a phase and we have proven the 
second assertion. 

We turn to the last claim. Let S and c be as defined in Eq. 
( [Tol l. Using the commutation relations Eq. (|4]i and the fact 
that conjugation by unitaries leaves the center x(t)l of the 
Weyl representation invariant, it is easy to see that 5* must be 
an isometry in the sense that [Sa, Sb] = [a,b]. To proceed, 
consider the following calculation. On the one hand 

Uw{a)w{b)U'' = Uw{a + b,2-'^[a,b])U^ (26) 
= w{S{a + b),2''^[a,b])c{a + b), 

while on the other hand, 

Uw{a)w{b)U^ = Uw{a)U'^Uw{b) (27) 
= w{Sa)w(Sb)c(a)c{b) 
^ w{Sa + Sb,2-^[Sa,Sb])c{a)c{b). 

Comparing the last lines of Eqs. ( |26] | and dZTl i one finds that S 
must be compatible with addition in Z^" meaning that 5(a + 
b) = Sa + Sb. Because TLd is cyclic the preceding property 
implies that S is also compatible with scalar multiplication: 

S{\a) = ^(a + • ■ • + a) = S{a) + ■■■ + S{a) = \S{a). 

Hence S is linear and therefore symplectic. Lastly, again us- 
ing lines ( |26] ) and ( l27T i. we have that c(a + 6) = c{a)c{b) and 
conclude that c is a character. By Lemma l24l there exists an 
ao e V such that c( • ) = x([ao, S ■ ]). Thus: 

w{aQ)^{S)w{a) ^{S)'^w{aQ)'^ — w{ao)'w{Sa)w{—ao) 

= x{[ao, Sa])w(Sa) 
= c{a)w{Sa). 

□ 

B. Axiomatic Characterization of the Wigner function 

The discussion in Section IIIDI should suggest that Defini- 
tion|5]yields 'the' natural analogue of the original continuous 
Wigner function. However, to bolster that claim with more 
objective arguments, we establish that - at least in prime di- 
mensions - the form is virtually determined by the property 
of Clifford covariance (Theorem|7]i. 



Theorem 23. (Uniqueness) Let d be an odd prime. Let 
Q,V,7i be as usual. Consider a mapping W' that fulfills the 
following axioms. 

L (Phase space) W' is a linear mapping sending opera- 
tors to functions on the phase space V. 

2. (Clifford covariance) W' is covariant under the action 
of the Clifford group, in the sense of Theorem^ 

Then Wp{p, q) = XiWp{p, q) + \2f0r two constants A1.2. // 
further, 

3. (Marginal probabilities) W' gives the correct marginal 
probabilities, as stated in Theorem \6\3\ 

then W'{p,q) = W{p,q). 

Proof Consider an alternative definition p W'p of a Wig- 
ner function. Linearity implies the existence of a set of op- 
erators {A'{v)} such that W'{v) = d"" tr(A'(w)p). W is 
covariant under the action of the Weyl operators if and only 
if A'{v) = w{v)A'{0)w{v)K So the only de gree of freedom 
left in the definition of W is the choice of ^'(0). Again, one 
must require A'{Sv) = ij.{S) A{v) Ijl{S) if Theorem |7] is to 
hold. In particular, because the origin is a fixed point of any 
linear operation, A'{0) must commute with all /i(S'). 

As a consequence, the old, unprimed Wigner function 
Wa'{q) of A'{0) stays fixed under any symplectic operation 
S. Since any two non-zero points of V can be mapped onto 
each other by a suitable symplectic matrix S, Wa'{v) must be 
constant on all such points. So there are only two parameters 
free to be chosen: T4^^'(o) (0) and Wa'{o) {v), v ^ 0. Clearly, 
the set of all operators that comply with these constraints is 
spanned by 1 and A{0): 

A'{0) = X1I + X2A{0). (28) 

The above decomposition imphes the first statement of the 
Theorem. 

As for the second claim, choose an a € V. The pro- 
jection operator \a){a\ is invariant under the action of Weyl 
operators of the form 'w{p, 0). Thus, due to Clifford covari- 
ance, the Wigner function VF'^^ must be p-shift invariant: 

^I'a) (P + p'tQ) = W^i'q) (Pj i)- We required Theorem |6]3] to 
hold, hence 

J2w(^^{p,0) = d" W'{0,0) = 6a,o- 

By Eq. (|28]l and Theorem iE it follows that W'{Q,0) = 
d-" {\i + A2(5a.o), yielding Ai = 0, A2 = 1. □ 

C. Characters and Complements 

Consider a space R = TL^ with a bilinear form ( • , • ) : 
R X R ^ TLd- For any s G R the function r 1-^ x((s, r)) de- 
fines a character of R. The form is said to be non-degenerate 
if (s, • ) ^ (s'j ■ ) for distinct ,s, s'. The two spaces we are 
concerned with are Q with the canonical scalar product and 
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V with the symplectic scalar product. Both can easily be 
checked to be non-degenerate. 

The following lemma states a basic fact about spaces with 
non-degenerate forms. We repeat it for completeness. 

Lemma 24. Let R = TH^ with non-degenerate bilinear form 
( • , • ). Any character of R is of the form C(r) — x([Si '']) 
for some unique s (z R. 

Proof. Addition gives V the structure of a finite abelian group. 
Therefore, y ~ y*, as is well-known (see e.g. Ref. [29]). So 
there are \V\ different characters of V, but equally many of 
the form x([w, •])• □ 

If d is prime and M a subspace of V, the well-known rela- 
tion dim AI + dim — dim V holds ['36']. It is, however, 
no longer true in the general case. A counter-example can be 
constructed along the same lines as in Section HTCl Still, an 
analogue exists as demonstrated below. 

Theorem 25. Let R = HI^ with non-degenerate bilinear form 
( • , • ). If M denotes a subspace of R, then the 'complemen- 
tarity relation' \M\ \M^\ = \R\ holds. 



Proof We will show that 



{VI MY 



(29) 



For m G M^, the relation [v] i— > v\) defines a char- 

acter of V/M, as can easily be verified. Let us denote the map 
TO 1-^ x{[m, •]) by ti. 

Conversely, given an element ^ of {V/M^ , v ^ C(b]) 
is a character of V. By Lemma |24] there exists a unique 
w ^ V such that C{[v]) = xii'^^'i'])- If G M, then 
({[m]) = C([0]) = 1 and hence w € M-^. Using the no- 
tions just introduced, we can define t2 : (V/M)* by 

It is simple to check that L2 ~ l^^. In particular, li is 
invertible and Eq. (|29] l follows. 

With the help of Lagrange's Theorem, we can compute 



|M^| = (V/M)* =\V/M\ = \V\/\M\, 
which concludes the proof. 



□ 



Corollary 26. Let V, Q be defined as usual. Let M be an 
isotropic subspace ofV and S be any subspace ofQ. 

L (Maximally isotropic spaces) M is equal to its symplec- 
tic complement M-^ if and only if\M\ = d". 

2. (Characters of subspaces) Any character C of S can be 
written as C(s) — x{qs) for a suitable q G Q. 

Proof. Claim [T] follows immediately from Theorem |25] and 
the fact that isotropic spaces are contained in their symplectic 
complement: AI C AI^. 

We turn to the second statement. In Lemma [l4l we have 
argued that the characters of S which are expressable as xiQs) 
stand in one-to-one correspondence to cosets in Q/S^. But 
\Q/S^\ = \S\ and hence all characters are of that form. □ 



D. A geometric note 

The proof of the Main Theorem makes use of the fact that 
for any vector w G F, there exists a symplectic operation S 
that sends u to a vector of the form {p,0). Indeed, if d is 
prime, any two vectors are similar, in the sense that they can be 
mapped onto each other by a symplectic matrix. Technically, 
this is a trivial incarnation of Witt's Lemma (see Ref. IstIi for 
a formulation that is applicable in our context). 

Once again the non-prime case poses additional difficulties. 
Recall that the order of a u e is the least positive A G 
such that \ v = 0. It is easy to see that the order of a vector 
is left invariant by the action of invertible linear mappings. 
If d is a composite number (i.e. not prime), then V = H?^ 
contains elements of different orders which cannot be related 
by a linear operation. However, one might conjecture that any 
two vectors of equal order are similar This is the content of 
the following lemma. Some concepts used in the proof can be 
found in Refs. LSHH. 



Lemma 27. (Similarity) Let V = K^". Let ai,a2 ^ V be two 
vectors with the same order Then there exists a symplectic 
matrix S such that Sai = 02. 

Proof. We can slightly weaken the assumptions made about 
V. All we require for this proof is that is a finite K^-module 
with non-degenerate symplectic form [ • , ■ ] . It need not be of 
the form Z^". 

Let V G V he a vector of order d. As v ^ x(b: ']) 
implements an isomorphism, V V*, ord(x([w, •])) — 
ord(v) — d. There hence exists aw such that [w, w] = A 
has order d. Any such number possesses a multiplicative in- 
verse A^^ modulo TLd and hence w' — A^^ fulfills [w, w'] ~ 1. 
Vectors satisfying such a relation are said to be hyperbolic 
couples. Denote their span {{v, w'}) as H. 

Set V := H^.^y Theorem|25] |y | = \H\\V'\. Further, it 
is easy to see that n if = {0} and hence V = HQ)V', 
where Q denotes the orthogonal direct sum. We claim that 
the symplectic inner product is non-degenerate on V' . Indeed, 
suppose there is a non-zero v' G V such that [u', w'] = for 
all w' G v. Then, by definition of V\ [h, w'] = for all 
h E H and therefore v' would be orthogonal on all vectors 
of V . Hence such a v' cannot exist by the non-degeneracy of 
[•,•]■ 

Note that V' fulfills the assumptions made about V at the 
beginning of the proof and has strictly smaller cardinality. 
Thus, we can induct on \ V\ to obtain a decomposition 



of V in terms of two-dimensional subspaces spanned by hy- 
perbolic couples {vi,w^}. We arrange these vectors as the 
columns of a matrix S = {vi, . . . ,Vn,w[, . . . ,iju'j^). The con- 
struction of the couples {vi, w'^} ensures that S is symplectic, 
as can easily be verified. 

Now let ai , 02 G be two vectors with maximal order 
By the preceding discussion, there exists symplectic matrices 
Si having Ui as their respective first column. Clearly, then 
<S'2<S'f ^fli = 02- 
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Lastly, suppose ord(ai) = k < d. It is easy to see that 
— ktti /d are elements of V with maximal order. Further, if 
S maps a'l to a'2, then also ai to 02. □ 

Corollary 28. (Transitive action) Let \Mi,vi), \M2,V2) be 
stabilizer states. If their respective associated isotropic sub- 
spaces Ml , M2 are spanned by vectors of maximal order, then 
there exists a Clifford operation relating these state vectors. 

Proof. Let {mf' , . . . , rrin' }, i = 1, 2 be bases of Mi and M2 
respectively. Assume that all vectors have maximal order. It 
is simple to adapt the previous proof for constructing a sym- 
plectic matrix S sending rnf^ to mf^. □ 

E. Some properties of the phase space point operators 

Lemma 29. (Properties of the phase space point operators) 
The phase space point operators fulfill the following relations 

A{a) = w{2a)A{0), 
A{a)A{h) = w{2a-2b), 
tr{A{u) A{v) A{w)) = xi[v,u] + [u^w] + [w,v]). 

Further, ifU permutes the phase space point operators un- 
der conjugation 

UA{v)U'^ ^ A{v') 
for all w e V, then U is Clifford. 



Proof. Clifford co variance (Theorem |7]l implies A{a) = 
w{a)A{Q)w{a)]. Using Theorem 16 15 1 it is easy to see that 
A(0)w(a)yl(0) ^ w{-a) and ^(O)^ = 1. Hence 

A{a) ^ w{a)A{0)w{-a)A{0) A{0) = w{2a)A{0) 

proving the first relation. The second one follows. 

For the proof of the third equation, we abbreviate A{0) as 
A. Then 

ti- {A{u)A{v)A{w)) 
= tr{w{2u)Aw{2v)Aw{2w)A) 
= tr{w{2u)w{-2v)w{2w)A^) 
~ xib^-j "v] + [u ~ v,w]) tr(w(2(u — v + w))A) 
~ u] + [u, w] + [w, v]) tr{A{u — v + w)). 

It has been noted in Theorem |6l6| that phase space point oper- 
ators have unit trace, which concludes the proof. 

Lastly, suppose the action of U permutes phase space point 
operators. For any a <E V,we have 

Uw{a)U^ = Uw{2 2-\a~0))U^ 

= UA{a)UU^A{0)U'' 

= A{a')A{0') 

= w(2(a'-0')) 

for suitable a', 0' G V. Hence U maps Weyl operators to Weyl 
operators and is thus Clifford by definition. □ 
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